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Galactic interstellar emission contributes substantially to Fermi LAT observations in the Galactic plane, where 
the majority of Supernova Remnants (SNRs) are located. We have developed a method to explore some 
systematic effects on SNRs' properties caused by interstellar emission modeling. We created eight alternative 
Galactic interstellar models by varying a few input parameters to GALPROP, namely the height of the cosmic 
ray propagation halo, cosmic ray source distribution in the Galaxy, and atomic hydrogen spin temperature. 
We have analyzed eight representative candidate SNRs chosen to encompass a range of Galactic locations, 
extensions, and spectral properties using the eight different interstellar emission models. The models were fitted 
to the LAT data with free independent normalization coefficients for the various components of the model along 
the line of sight in each region of interest. We will discuss the results and compare them with those obtained 
with the official LAT interstellar emission model. 



1. INTRODUCTION 

Galactic interstellar 7-ray emission is produced 
through interactions of high-energy cosmic ray 
(CR) hadrons and leptons with interstellar gas 
via nucleon-nucleon inelastic collisions and electron 
Bremsstrahlung, and with low-energy radiation fields, 
via inverse Compton (IC) scattering. Such interstellar 
emission accounts for more than 60% of the photons 
detected by the Fermi Large Area Telescope (LAT) 
and is particularly bright toward the Galactic disk. 

In this paper, we present our ongoing effort to ex- 
plore the systematic uncertainties due to the model- 
ing of Galactic interstellar emission in the analysis of 
Fermi LAT sources, with particular emphasis on its 
application to the 1*^ Fermi LAT Supernova Remnant 
(SNR) Catalog. We compare the results of analyz- 
ing sources with eight alternative interstellar emission 
models (lEMs), described in Section [2l to the source 
parameters obtained with the standard model in Sec- 
tion [3l In Section m we discuss the future application 
of this method to the SNR Catalog. 



2. INTERSTELLAR EMISSION MODELS 

In order to estimate the systematic uncertainty in- 
herent in the choice of standard interstellar emission 
model (lEM) in analyzing a source, we have developed 
eight alternative lEMs. By comparing the results of 
the source analysis using these eight alternative mod- 
els to the standard model, we can approximate the 
systematic uncertainty therefrom. 



2.1. The Standard lEM 



The standard lEM for Fermi LAT data analysis 
was developed by the collaboration using the sim- 
ple assumption that energetic CRs uniformly pene- 
trate all gas phases in the interstellar medium. Un- 
der this assumption, the Galactic interstellar 7-ray 
intensities can be modeled as a linear combination of 
gas column densities and an inverse Compton (IC) 
intensity map as a function of energy. The gas col- 
umn densities are determined from emission lines of 
atomic hydrogen (Hi)^ and CO, the latter a surrogate 
tracer of molecular hydrogen, and from dust optical 
depth maps used to account for gas not traced by 
the lines. To account for a possible large scale gra- 
dient of CR densities, the gas column density maps 
were split into 6 Galactocentric rings using the emis- 
sion lines' Doppler shift. The IC map is obtained us- 
ing GALPROP^ to reproduce the direct CR measure- 
ments with a realistic model of the Gala ctic interstel- 
lar ra diation field (ISRF), as was done in lPorter et al\ 
J2QQ8( . To account fo r some extended remaini ng resid- 
uals, notably Loop I Casandiian et al. [2009] and the 
so-called Fermi bubbles Su et al. [2010], the standard 
lEM includes them as additional templates. 



■"^H I column densities are extracted from the radio data using 
a uniform value for the spin temperature (200 K). 

^The GALPROP code has been developed over se veral 
years, starting with, e .g. Mo skalenko and Strong] |1998[ and 
IStrong and Moskale"nka ^1998i] . 
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This lEM, along!: w ith sources in the 2FGL Cata- 
log 'Nolan et al\ |2Q12 ] and an isotropic intensity ac- 
counting for the extragalactic 7-ray and instrumental 
backgrounds, was fit to two years of LAT data. This 
yielded best fit values of the linear combination co- 
efficients which can be interpreted as gas emissivities 
in the various Galactocentric rings and a renormaliza- 
tion for the IC model as a function of energy, as well 
as the spectrum of the isotropic component. The ra- 
tio of CO to twice the H i emissivities and dark gas to 
Hi emissivities are proportional to the CO-to-H2 and 
dust-to-gas ratios, respectively, under our simple as- 
sumption. Another explanation of this method of de- 
composing the 7-ray sky to create the stan dard model 
may be found in, e.g. iKatagiri et al\ |2Qll( . The stan- 
dard lEM is distributed as a cube summed over the 
components which predicts the intensities of Galactic 
interstellar 7-ray emission in a grid of directions and 
energies with its accompanying isotropic model. Fur- 
ther description and details are available at the Fermi 
Science Support Center^. 



2.2. Alternative lEMs 

To explore the uncertainties related to the modeling 
of interstellar emission we generated eight alternative 
lEMs, in particular probing key sources of systematic 
uncertainties by: 

• adopting a different model building strategy 
from the standard lEM, resulting in different 
gas emissivities, or equivalently CO-to-H2 and 
dust-to-gas ratios, and including a different ap- 
proach for dealing with the remaining extended 
residuals; 

• varying a few important input parameters for 
building the alternative lEMs: atomic hydrogen 
spin temperature (150 K and optically thin), CR 
source distribution (SNRs and pulsars), and CR 
propagation halo heights (4 kpc and 10 kpc); 

• and allowing more freedom in the fit by sepa- 
rately scaling the inverse Compton emission and 
H I and CO emission in 4 Galactocentric rings. 



The work in lAckermann et al\ |2012a( , using the 
GALPROP CR propagation and interaction code, was 
used as a starting point for our model building strat- 
egy. The GALPROP output intensity maps associ- 
ated with Hi, CO, and IC are then fit simultaneously 
with an isotropic component and 2FGL sources to 
2 years of Fermi LAT data in order to minimize bias in 
the a priori assumptions on the CR injection spectra 



and the proton CR source distribution. The inten- 
sity maps associated with gas were binned into four 
Galactocentric annuli (0 — 4 kpc, 4 — 8 kpc, 8 — 10 kpc 
and 10 — 30 kpc). The spectra of all intensity maps 
were individually fit with log parabolas to the data, 
to allow for possible CR spectral variations between 
the annuli for all Hi and CO maps while the IC fit 
accounts for spectral variations in the electron dis- 
tribution. We also included in th e fit an isotropic 
template and templates for Lo op I ICasandiian et al\ 
[200a] and the Fermi bubbles ISu et al\ [2Q10J . The 
te mplate for Loop I is based on the geometrical model 
of I Wollebe n [2007] while bubbles are assumed to be 
uniform with edges defined in spherical coordinates 
hy R = Rq\ sin^l, where is the polar angle. 



^http://fermi. gsfc.nasa.gov/ssc/data/access/lat/ 
ModeLdetails/Pass7_galactic.html 



Ackermann et al\ [2012a| explored some systematic 
uncertainties by varying input parameters. The Hi 
spin temperature, CR source distribution, and CR 
propagation halo height were found to be among those 
parameters which have the largest impact on the 7- 
ray intensity. The values adopted in this study to 
generate the eight alternative lEMs were chosen to be 
reasonably extreme; we note that they do not reflect 
the full uncertainty in the input parameters. Sepa- 
rately scaling the Hi and CO emission in rings and 
the IC emission permits the alternative lEMs to bet- 
ter adapt to local structure when analyzing particular 
source regions. Figure [1] shows the relative difference 
between the standard model and one of the alterna- 
tive models (Lorimer CR source distribution with a 
4 kpc halo height, and 150 K Hi spin temperature). 
Differences are particularly large along the Galactic 
plane, where SNRs are located. 

Finally, we note that this strategy for estimat- 
ing systematic uncertainty from interstellar emission 
modeling does not represent the complete range of sys- 
tematics involved. In particular, we have tested only 
one alternative method for building the lEM, and the 
input parameters do not encompass their full uncer- 
tainties. Further, as the alternative method differs 
from that used to create the standard lEM, the result- 
ing uncertainties will not bracket the results using the 
standard model. The estimated uncertainty does not 
contain other possibly important sources of system- 
atic error, including uncertainties in the ISRF model, 
simplifications to Galaxy's geometry, small scale non- 
uniformities in the CO-to-H2 and dust-to-gas ratios 
and Hi spin temperature non-uniformities, and un- 
derlying uncertainties in the input gas and dust maps. 
While the resulting uncertainty should be considered 
a limited estimate of the systematic uncertainty due 
to interstellar emission modeling, rather than a full 
determination, it is critical for interpreting the data, 
and this work represents our most complete and sys- 
tematic effort to date. 
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(PL) spectral mode l either replaces th e closest non- 
pulsar 2FGL source iNolan et aU 2Q12J within the ra- 
dio size or is positioned as a new sou rce at the loc ation 
determined from radio observations iGreenI |2QQ9( . All 
other 2FGL sources within the radio size which are 
not pulsars are removed from the source model. We 
fit the centroid and extension of the SNR candidate 
disk as well as the normalization and PL index for 
the source of interest and the five closest background 
sources within 5° with a significance of > 4 a. This 
procedure balances the number of degrees of freedom 
with convergence and computation time requirements. 



Figure 1: The position of the eight candidate SNRs 
used in this analysis are overlaid on a map of relative 
difference between the standard lEM and one of the 
alternative models. The alternative model selected 
for this image has a Lorimer source distribution, a 
halo height of 4kpc and a spin temperature of 150 K. 
We plot the difference between the models' predicted 
counts divided by the square of the sum of the 
predicted counts so the map is in units of sigma. 
The hardness of the SNRs' spectra is in two 
categories: hard (purple) and soft (black). SNRs 
detected as extended with Fermi are shown as circles 
while point-like are shown as crosses. 



3. ESTIMATING lEM SYSTEMATICS 



3.1. Analysis Method 

We developed this method for estimating the sys- 
tematics from the interstellar emission model using 
eight candidate SNRs chosen to represent the range of 
spectral and spatial SNR characteristics in high and 
low lEM intensity regions. Figure [1] shows the can- 
didate SNRs' location on the sky, illustrating their 
range of Galactic longitude. The color indicates those 
candidates with a hard or soft index and the shape 
the extension (pointlike or extended). The SNR can- 
didates are overlaid on a map of the relative difference 
between the standard lEM and one of the alternative 
lEMs described in Section [2j 

We use the same analysis strategy to obtain all SNR 
candidates' Fermi LAT parameter values with both 
the standard and all eigh t alternative lEMs on 3 ye ars 
of P7_V6SOURCE data lAckermann et al\ J2Q12bj in 
the energy range 1 — 100 GeV. We applied the standard 
binned likelihood method^, treating sources as follows. 
For each of the eight candidate SNRs, an extended 
source initially of the radio size and with a power law 



L=-30' 



L=-60* 




^The stand ard Fermi LAT analysis description and tools can 
be found here: 'littp://fermi. gsfc.nasa.gov/ssc/data/analysis/ . 



Figure 2: Schematic representation of the Hi and 
CO rings used for the split alternative lEMs crossed 
by lines of sight at various Galactic longitudes. The 
ring boundaries lie at 0, 4, 8, 10, 30kpc from the 
Galactic center while the red dot marks the sun 
position at 8.5 kpc. The figure is not to scale. 



To generate results for the source of interest with 
each diffuse model, we fit the sources' model to the 
data with either the standard model or one of the 
eight alternative lEMs. In the case of the standard 
Fermi LAT lEM, we allow the normalization to vary 
and fix the accompanying standard isotropic model's 
normalization. For each of the eight alternative mod- 
els, we use the corresponding isotropic model fixed to 
its value resulting from the fit to the all-sky data (see 
Section [2]). To better understand the effect of allowing 
freedom in the H i and CO rings, we fit the alternative 
models in two ways: either with the rings' normaliza- 
tions free ("split" models) or with the rings summed 
together, as given by the all-sky fit (see Section[2j), and 
only the total normalization free ("summed" mod- 
els). The summed alternative lEMs are thus closer 
to the standard lEM. For the split alternative lEMs, 
as shown in Figure [21 not all rings are crossed by all 
lines of sight. We thus fit only the two innermost Hi 
and CO rings crossed by the fine of sight to our region 
of interest. The IC template is also free to vary while 
the isotropic component remains fixed. 
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(a) Flux for the eight candidate SNRs' from 1 - 100 GeV. (b) Index for the eight SNR candidates from 1 - 100 GeV. 

Figure 3: Results for each candidate SNR, averaging over the eight alternative lEMs separately for split (red) 
and summed (green) component models compared to the standard model solution (black). The error bars for 
results using the alternative lEMs show the maximal range of the values given by the 1 a statistical errors. 



3.2. Results for SNRs' lEM Systematics 

To compare the results obtained using the eight al- 
ternative lEMs with the standard model results, we 
average each parameter's eight values from the alter- 
native lEMs. Figure [3] shows the values for the flux 
and index from fitting the data with the alternative 
lEMs with the rings either split or summed. These are 
then plotted along with the standard model results for 
all eight SNR candidates studied. We conservatively 
represent the allowed parameter range with error bars 
showing the maximal range for the alternative lEMs 
1 a statistical errors. 

Figure [3] shows that the variation in value of the best 
fit parameters obtained with the alternative lEMs is 
larger than the 1 a statistical uncertainty. The im- 
pact of changing the lEM on the source's parame- 
ters depends strongly on the source's properties and 
location. As expected, the parameter values for the 
source of interest are generally closer to the standard 
model results for the alternative lEMs with compo- 
nents summed rather than split. In many cases, the 
allowed parameter range represented by the 1 a statis- 
tical errors for each of the alternative lEMs is larger 
with the components split than summed. Also as 
noted earlier, the alternative lEM results do not as a 
rule bracket the standard model solution. We observe 
that some of the largest differences between the stan- 
dard and alternative lEM results for a single source 
are frequently associated with sources coincident with 
templates accounting for remaining residual emission 
in the standard lEM (Section [2]). 

SNR G347.3-0.5 proves an interesting source for un- 
derstanding the impact nearby source (s) can have on 



this type of analysis. In particular, our automated 
analysis finds a softer index and a much larger flux 
for SNR G347.3-0.5 than th at obtained in a dedicated 
analysis lAbdo et all |201l( . Since the best fit radius 
(0.8°) is larger than that the X-ray data indicates 
(0.55°), the automated analysis's disk encompasses 
nearby sources that are only used in the lAbdo etal\ 
|2Qll( model. Including this additional emission also 
affects the spectrum, making it softer in this case than 
that found in the dedicated analysis. Given Fermi 
LAT's both increasing point spread function and num- 
ber of sources with decreasing energy as well as the 
predominance of diffuse emission at lower energies, we 
note that nearby sources may play a greater role if 
extending this method below the 1 GeV minimum en- 
ergy examined here. 



3.3. lEM Input Parameter Comparison 

To identify which, if any, of the three lEM in- 
put parameters (H i spin temperature, CR source dis- 
tribution, and CR propagation halo height) has the 
largest impact on the fitted source parameters, we 
marginalize over the other parameters and examine 
the relative ratio of the averaged input parameter 
values to the values' dispersion. For a fitted source 
parameter a, such as flux and a GALPROP input 
parameter set P = {^,j}, e.g. spin temperature 
Ts = {150K, 10^ K}, this becomes: 



< ai> 



max{aa,ijO-a,^ 



(1) 
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Figure 4: Schematic representation of the method for 
determining the impact of a given GALPROP input 
parameter (e.g. spin temperature) on a source 
parameter (e.g. flux) by marginahzing over the other 
input parameters. We take as the figure of merit the 
ratio of the difference in averages for 
Ts = {150 K, 10^ K} divided by the maximum RMS. 



where aa is the rms of the parameter a for a given in- 
put parameter value P. Figure 3] shows this schemat- 
ically. A ratio > 1 implies that changing the selected 
input parameter has a greater effect on the fiux than 
all combinations of the other input parameters. 

In Figure [5] we plot this ratio for each of the alter- 
native lEM's input parameters for each of the eight 
SNR candidates, along with the average over the SNR 
candidates, separately for the split and summed com- 
ponents. While the spin temperature has the largest 
effect for the split alternative lEMs, the CR source 
distribution also becomes relevant with the summed 
alternative models. In light of this and as none of the 
parameters shows a ratio significantly greater than 1 
for all the sources tested, we conclude that none of 
the input parameters has a sufficiently large impact 
on the fitted source parameter to justify neglecting 
the others. 



4. FUTURE APPLICATIONS 

In this work we explored the effect of using alter- 
native interstellar emission models on the analysis of 
LAT sources. As the Galactic interstellar emission 
contributes substantially to Fermi LAT observations 
in the Galactic plane, the choice of lEM can have a 
significant impact on the parameters determined for 
a given source of interest, as demonstrated with eight 
SNR candidates. To estimate the reported error we 
currently use only the most conservative extreme vari- 
ation of the source of interest's output parameters. 
We are finalizing our definition of the systematic er- 



ror using this method, including through comparison 
of the present estimate with previous methods' esti- 
mates, typically found by varying the standard lEM's 
normalization by a fraction estimated from neighbor- 
ing regions. Although our current method represents 
the uncertainty due to a limited range of lEMs, it 
plays a critical role in interpreting the data and rep- 
resents the most complete and systematic attempt at 
quantifying the systematic error due to the choice of 
lEM to date. 

As the majority of SNRs lie in the Galactic plane, 
coincident with the majority of the Galactic interstel- 
lar emission, this method is particularly pertinent to 
analyses such as that underway for the 1*^ Fermi LAT 
SNR Catalog. Figure [3] shows that the fiux and in- 
dex can vary greatly for our eight representative SNR 
candidates, depending on the source and local back- 
ground's specific characteristics. Given these differ- 
ences, we plan to use this method to estimate the sys- 
tematic uncertainty associated with the choice of lEM 
on the full set of SNR candidates in the catalog. Such 
error estimates will allow us to, among other things, 
more accurately determine underlying source charac- 
teristics such as the inferred composition (leptonic or 
hadronic) and particle spectrum. 

Other classes of objects such as pulsar wind nebulae 
and binary star systems also lie primarily in the plane 
and are likely to be strongly affected by the choice 
of lEM. We are thus generalizing this method in or- 
der to be able to apply it to the study of Galactic 
plane sources generally. Another possible extension 
to this method is extending it to energies < 1 GeV, 
where the interplay between the Galactic interstellar 
emission model and background sources must be care- 
fully examined. By more faithfully accounting for the 
systematic uncertainty of our model components we 
will be better equipped to draw less biased conclu- 
sions from our data. 
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Figure 5: The impact on the candidate SNRs' flux of each of the alternative lEM input parameters, 
marginalized over the other GALPROP input parameters, is shown relative to the figure of merit for the other 
input parameters (source distribution, halo height, and spin temperature). We calculate the figure of merit 
(Eq[T]) separately for the alternate lEM components fit separately (left) and summed (right). The large open 
cross represents the average figure of merit over all SNR candidates. As no alternative lEM input parameter 
has a figure of merit significantly larger than 1, no input parameter dominates the fitted source parameter 
sufficiently to justify neglecting the others. 
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